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EFFECT OF RESONANCE -OB LATENESS COUPLING 
ON A SATELLITE ORBIT 
C. A. Wagner 

Goddard Space Flight Center 
Greenbelt, Maryland 

ABSTRACT 

Second order effects of the coupling between geopotential resonance and 
oblateness on a satellite orbit are calculated. These effects arise from the 
interaction of resonance with the secular changes of the orbit's node, perigee 
and mean anomaly. They have the same period and phase as first order 
resonance perturbations. But their amplitudes are proportional to the square 
of the period knd dominate the first order effects as the orbit becomes com- 
mensurate. A striking example of this coupling is seen in the 18 day resonance 
variation of the node of the orbit of the first earth resources technology satellite. 
Analysis of this one arc second (31m) variation 5 delds a strong 14th order 
constraint to the geopotential field for odd degree terms. This constraint 
(lumped coefficient) is poorly predicted by current models. 
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INTRODUCTION 

The analytic perturbations of a Kepler orbit by the geopoetntiai have been 
extensively discussed in the literature [i.e.; Kaula, 1966] . The zeroth order 
perturbations are the integration of the usual non-Unear Lagrange planetary 
equations: 

E = f (E, J, (1) 

with constant right hand side (except for the steady progression of the mean 
anomaly). In Equation (1), E is a vector of six Kepler elements. The and 
So are the usual geopotential coefficients (but fully normalized; Kaula (1966), p, 7) 

'V, rn 

defined from the representation of the potential as: 


V =it 
* r 


00 /C’ 

-EL 

a 2 m “ 0 


(a/r)'^ P^^(sin0) [C^ 


cos m\ + So sin m\] 

III 'vm 


where r,<p , k are the satellite's distance from the earth's center of mass, 
geocentric latitude and longitude; fi is the Gaussian gravitational constant, a^ 
is the earth's mean equatorial radius and the ^ are the usual (but fully 

normalized) associated Legendre polynomials. At this initial level of approxi- 

0 

mation the well known secular effects on the node (H), perigee (w ) and mean 
anomaly (M) due to oblateness (C^ ^ ) dominate since this term is about one 
thousand times larger than the others. These zeroth order secular perturba- 
tions are: 


1 
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3/x^'^^ Q • /S cos I 

2(1 - e^)^ 


3^ 


l/2f 


■' 2,0 


, /g ,,2 q -7/2^1 _ C ^,-,^2 


a 

e 


!(1 - 5 cos'' I) 


( 2 ) 


"o = 


4(1 - 6“') 


2^2 


\ = 1 ^ [3 cos2 I - 1] , 

4(1 - e2)3/2 


where a, e and I are the orbit’s semimajor axis, eccentricity and inclination. 
Zeroth order perturbations of short period (less than an orbit revolution) also 
exist, again dominated by the Cj o term. However, these are usually calculated 
with respect to the secularly processing ellipse defined by Equation (2) and 
are thus first order perturbations. 

First order variations; integrations of Equation (1) incorporating the zeroth 
order effects (of the precessing ellipse) in their right hand sides, are also well 
known. The best example is the "odd" zonal (m = 0) oscillation of long period in 
e, i,<iJ , O and M with frequency w . 

In general, the first order perturbations have frequencies: 

ij) = (^ - 2p) «{| + (-t - 2p + q) M + tn(ng + 9^), 

where p and q are additional indices related to the representation of the geopo- 
tential in terms of the Kepler elements, 6^ is the rotation rate of the earth and 
M = n + Mq , vv^ere n = g[-3/2 , jn particular, longitude harmonic (m / 0) 

resonance occurs when’^= 0. 
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As is well known, the first order perturbations are of order "" 

(C|_^ , But of course, first order resonance perturbations are not 

adequate as>/(-0 (i.e., Gedeon, 1969). Evenfor so called ’’shallow resonance," 
where the character of the motion is still sinusoidal, a second order perturbation 
in the mean anomaly, of order , dominates the first order variation. 

This particular second order perturbation arises from the interaction 
(coupling) of the first order resonance variation of the semimajor axis with the 
two body mean motion. Following Kaula (1966), the element rates are expanded 
in a Taylor series with respect to the element perturbations: 




For the semimajor axis, the first order variation, of order J^^(n/'/'), for a 

< 

particular term ('t, m, p, q) is: 

Aj 3 = (n/i^) a(aya)^ 2F^^^(I) - 2p + q) 

where F and G are inclination and eccentricity functions and; 

n't “ m even 


(5) 


( 6 ) 


tn p q 




/mpq 


1 

eP 
/ 1 

even 

cos Ip + 



odd 




s in 0. 


I m odd 

[The characteristic longitude is i/; = (•i - 2p)co + (-^ - 2p + q) M + m(G - )]. 

Writing Equation (5) for the variation of the mean anomaly; 
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M = n + :^A, a, (7) 

oa 

the second term on the right gives the interaction with the linear perturbation 
of Equation (6): 

AjM = (- 3n/2a) Aj a. 

The integration of this "second order" effect with respect to [letting dt = 
yields: 

AjM = - 3(a^/a/ TG(l - 2p + q) S - (AM)^ S (8) 

where 

S = / Sd^. 

The second order (resonance) interaction perturbations to be discussed 
in this paper are entirely analogous to A M except that they are of order 
^ 2,0 * Though 10"^; \p << 1 in resonance and these effects 

will dominate the linear perturbations of order (n/^) an ^ - 0. 

ANALYSIS 

i 

The resonance perturbations of the node, perigee and mean anomaly are 
written as a sum of direct (linear) and indirect effects: 


A(n, fit), M) — Aj (fl, Ci>, 


BHo’ ^o> ^ . 

^ ~ A| 3 , c, I 

oa, e, I 


+ AjMCTWO BODY INTERACTION) 


(9) 
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Using the zeroth order secular rates of Equation (2); 


30 


3a 


" = - 2^1/2 c, 


2,0 


/5 COS I a^a“^^^/4(l - e^)^, 


— -^= 6/Li'^^C, Q * 1^5 e cos la^ (1 - e^)^, 

de 

30j, - 3ij}^^C2 o • v^sin I 

2(1 - e 2)2 


— = - 21fj}^^C2 0 • 1/^(1 - 5 cos2 I) a2a'®/V8(l - e^)^, 
3a ’ 

^ 0 • A e(l - 5 cos2 I) a^ /(i _ e2)3, 

00 ’ 


3<i„ 

31 


= Q • v^sin I cos I a^a"’^^/2(l - e^)^. 


^ = 21m‘"'“Cj.„\- (3 cos" I - 1) aj a-«/"/8(l - e")"/", 

O 3 


^ > V/2 c, • /^e(3 cos2 I . 1) a2 a-®/=^./4(l - , 

3e -ir 


and 


315^ 

hl- 


= Cj 0 • sin I cos I a^a''’^^/2(l - e^) 


2\3/2 


The first order (resonance) perturbation Aj^ a has already been given in 
Equation (6). The other direct resonance variations are: 
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Aj 6 = Aja{(l - [(1 - - (^ - 2p)]}/2ae 

Aj I = Aj^a {('i - 2p) cos I - m} /2a(l - s in I (11) 


Recall that resonance occurs when M = m/k, or a rational fraction of the 
earth's rotation rate. Here, k = (-t - 2p + q) = 1, 2, 3, . . . In these circum- 
stances, 0, A resonant family is characterized by m and k having all the q 
indices as species (or members). A specification of m and q however is sufficient 
to denote a particular species (or member) when the orbit's mean motion is 
known. It is noted that the strongest resonances usually occur in the families 
for k = 1 when q = 0 since the G functions are of order e ^ L 

To show the relative dominance of the direct and interactive terms, recall 
[from Kaula, 1966 p. 40l [that the direct effect in the node is given as: 


where 


AjQ = 


(n/^) (a/a) dF/dl GS 
(1 - sin I 


= (AD)i S, 


S = J sm d^ 


(12) 


The indirect, resonance -oblateness coupled perturbation is evaluated from 
Equations (6), (9), (10), and (11). For the integral, the change of independent 
variable d\p =^dt is made and all parameters except are assumed constant. 
Then 
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S(n/0)q,o 

(1 - (a/a^)2 


{2(n/0) FGk/(a/a^)^} 


- 2 1 cos I 


+ L£^H [(1 - ,e2)l/2 [(1 _ g2)l/2 _ (k - q)]] - 3 [(k - q) cos I - M\ ] _ 
(I - e2) 4(1 - e2)i/2 J 


Combinirg (12) and (13); 


An (TOTAL) = S GF J (^F/3I)/F 

(a./aj^ \(l-e2)l/2 sin I 


2(n/0) Cj 0 • /S k 

(l-e2)2 (a/aj2 


- 21 cos I 
4 


3 cos I _ ^ 2 y /2 
(1- e2) 


[(1-02)1/2 


- (k - q)] 


3 [(k - q) cos I - m] 

^ (1 - e2)l/2 


(14) 


Consider a resonant member of the family k = 1, m = 15 (M = 15 revolutions/day, 
or n ~ 5300 °/day). The dominant resonant effects are generally due to the 
lower degree terms where - m is small. Their sinusoidal inclination 
functions have characteristic wavelengths of 2^/-^ - m + 1 (Allan, 1973). 

T 3 T)ically then, in about one radian the chan g e in F is equal to F, or 3F/3I = 

0(F). Better, for a sinusoidal function F with wavelength of ^ , the value of 
3F/3 I (rms) is 2 tt/X times the value of F (rms). (A more exact "average" is 
derived in Appendix A). Thus for a close satellite orbit the order of the direct 
term will be about 4 (taking sin I into consideration) while that of the indirect 
term will be (Wi^) x lO ^ x lo or 55“^*^/ In other words, at a resonant 
period of about 25 days, the two effects are about the same for the average 
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close orbit. This ’’shallow resonance" period is somewhat longer than the majority 
of the useful geodetic orbits. The typical geodetic orbit has a dominant resonance 
period of about 4 days. But even with indirect effects 1/6 of the direct, these 
interaction terms cannot be ignored. 

Returning to the calculation of the second order oblateness-resonance effects 
in terms of the first order perturbations (A^ajA^ e, A^I) from Equations (9), 

(10) and (11); these are combined with the direct effects to give the total 
resonant variations as; 


An = /(An)j + 


(n4)c^,’ys 

(a/a )2 (1 - e2)2 


-21 


cos I(Aa)j^/a 


6e cos I(Ae), -i 

+ 1 -i sin I(AI), 

<1 - ^ 


Ao) = + 


(n/ipyc^ Q ■ v/S 

(1 - e 2)2 


-21 


(1-5 cos^ I) (Aa)j/a 


3e(l-5cos2l) 15 sin I cos I 

+ + (Al)j 

(1 - e2) 2 


S, 


and 


Am = (am)i + (am)j + 


(n/^) 

(a/aj2 (1 ^ e2)3/2 


21(3 cos^ I - 1) (Aa)j/a 


9e(3 cos^ I ' 1) 9 sin I cos I(Al)j 

^ - 


4(1 - e^) 


S, 


(15) 
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where: 

(Aw)j = AjO)/S, (AM)j = \U/S, (Aa)j ^A^ a/S 

(Ae)i =Aje/S, (Al)j =Ajl/S, 

and the Aj variations are the linear resonance effects given by Kaula, 1966, p. 40, 

An important point to note from Equation (15) is that both the phase and 
the frequency of the coupling effects are identical to those of first order. They 
will therefore be ordinarily indistinguishable from the linear resonance pertur- 
bations on a single orbit. But incorrect resonant coefficients will result in 
gravitational determinations from many orbits which have significant coupling 
and employ only Linear effects. Generally, the error in the mean anomaly and 
argument of perigee is not serious. This is because the in-plane resonance 
variation is usually dominated by the stronger second order resonance coupling 
with the two body mean motion which have always been included in analyses 
(i.e., Gaposchkin and Lambeck, 1970). However, the resonance out-of-plane 
from the variation of the node may have serious errors if only linear analysis 
is used as in Gaposchkin and Lambeck, 1970 and Gaposchkin, 1973. Table 1 
presents these maximum linear and total variations of the node for the 21 satellite 
orbits (excluding ERTS) used in the Smithsonian Standaixi Earth's (see also 
Gaposchkin, 1966). Almost half show more than 10% change when the indirect 
effects are included. 
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DISCUSSION AND RESULTS 


The analysis here is a straightforward application of perturbation theory. 

In fact, Kaula 1966, p. 49 has worked out two such second order variation; one is 
a long period variation of the node from odd zonal-Jj interaction. These effects 
are entirely analogous to those here. Born (1974) has worked out coupling effects 
of order n/^ f from the interaction of periodic resonant terms with 
themselves. This was done for the analytic theoi^y of the Mariner 9 Mars orbiter 
(n 2 revs/day) which had a shallow resonant period of 18 days. For the close 
earth satellite these resonant-resonant interactions are of the order of 10 ^ 
less than the resonant-oblateness ones. For Mars, the longitude gravity field 
is more than an order of magnitude stronger than on earth and the high Mariner 9 
orbit further diminishes the obleteness effects. As a result, the resonant- 
resonant effects are roughly comparable to the resonant-oblateness ones. More 
surprising, even though the Mars oblateness is almost twice earth’s, the indirect 
effect on the Mariner 9 node is only about 2% of the direct. This is due to the 
height of the orbit and a significant cancellation of indirect terms. 

For example, I use the following parameters for the Mariner 9 orbit about 
Mars: 

a = 12631 km(~2 rev’s/day), e = 0.62, I = 64.8", 360" /<^ = 18.3 days 
with respect to the term (2, 2, 1,1) 10® C 2 ,o^“ " 1960//5, 10® ^ ” 96.7. Then, 

Equation (12) gives the amplitude of the direct perturbation of the node as; 

lAjDl = .02580". 

This agrees well with Born (1974). On the other hand the indirect effect from 
Equation (13) is calculated as: 

10 
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lAjfil = .00051“ (17) 

This amplitude is somewhat greater than the amplitude of the periodic part of 
the resonant-resonant interaction (Born, 1974). It also appears to be sufficient 
to explain the remaining discrepancy of the Mariner 9 theory (for the node) with 
numerical integration. 

ERTS RESULT FOR VARIATION OF THE NODE 

The ERTS 1 spacecraft was intensively tracked by Unified S-Band range 
and range rate stations (worldwide) in 1972 and 1973. From this tracking 
history, precise, almost daily osculating elements have been obtained (Arthur 
Fuchs, Goddard %>ace Flight Center). These have been converted to mean 
elements by numerically averaging minute Keplerian ephemerides over each 
day of tracking, with respect to a processing ellipse. Before averaging, the 
osculatii^ vectors were first reduced of most short and medium period terms 
analytically (Douglas; Marsh and Mullins, 1972). 

Figure 1 (circled points) shows a residual variation of the mean node in a 
50 day period in late 1972. These residuals are the "observed” mean nodes minus 
calculated values from a trajectory computed with all the long term effects of 
the zonal geopotential, the sun and moon (including earth tides)., radiation pressure 
and atmospheric drag (Wagner, Doi^las and Williamson, 1974). No resonant 
geopotential effects were Included on this trajectory. 
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The ERTS orbit’s ground trace repeats in 18 days (by design) which is seen 
to be the period of this variation. The amplitude is only about 30 m but well 
defined. The ERTS orbit parameters are; n - 14 rev’s day, e = .001 and I = 99.1°. 
It is natural, therefore, to ascribe this variation to a resonance with geopotential 
terms of order 14. In particular, from Equation (3), 360° /p = 18.0 days, when 
= (o> + M) + 14(n +0^). This characteristic rate holds for all odd degree resonant 
terms (-t , m, p, q) of order 14 wdiere q = 0. By almost two orders of magnitude 
this f ami ly has the strongest effect on the node. (See Table 2.) 

C. A. Wagner and S. M. Klosko, unpublished, have analysed the variation 
of the semi-major axis, inclination, node and eccentricity of the ERTS orbit 
by the same semi-numerical trajectory program which calculated the residuals 
in Figure 1. In this analysis only a pair of two odd and even degree coefficients 
of 14th order was used, with a single 28th order pair. They found the pair of 
two ’’lumped" odd degree, 14th order coefficients which best fit the mean element 
data to be: 


10® (C, S)jg = - 23.1 ± 0.7, - 97.5 ± 0.7 


(18) 


10® (C, = 37.0 ± 0.8, 138.4 ± 0.8 

Using these values in Equation (12) (with the sum of the two terms taken) or 
scaling to the values in Table 2, the analytic variation of the node considering 
only first order effects is: 


Af2" = - 0.347 sinp - 0. 149c os i/; 


(19) 


12 



15 


The corresponding variation including the 2nd order effects of oblateness-resonant 
interaction is (from Equation (15)): 

AD" = - 0.696 s in i/j - 0.238 cos 

The longitude 4< in the 50-day period shown in Figure 1 is nearly linear in time: 

v/;° = - 11.7 + 18(MJD - 41622.33). (21) 

The two variations are plotted in Figvire 1. Clearly, the second order varia- 
tion is more faithful to the data. 

However, it should be apparent that any fixed biases in the trackers such as 
station coordinate and constant calibration errors can simulate this resonant 
effect (i.e., Allan, 1973, p. 223). Such tracking depends on the orbit-station 
configuration and would repeat in the same period as the ground track, which is 
the resonant period. It is difficult to determine the likelihood or m^nitude of 
such an effect over a complex of tracking stations and operating times. How- 
ever, with regard to station position error (center of mass coordinates), the 
ERTS trackers should be known to better than 10m individually. Any residual 
bias in the node should be considerably less than this in orbits determined from 
all the stations. It should be added that such a bias would affect all geodetic 
results from resonant satellites, but especially those with effects less than 10m. 

A considerable amount of the shallow resonant information in current satellite 
determined fields is only available at a level below 10m (i.e. Table 1 and 
Lerch et al., 1974). 
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Finally.it is instructive to write out the full resonant constraint for the nodal 
variation of the ERTS orbit as defined by Equation (20). The lumped coefficients 
of this constraint are merely the coefficients of the sine and cosine terms in 0 ex 
panded as a sum of all the relevant 14th order harmonics. The ’’observed" 
lumped coefficients can then be compared to calculated values from more com- 
plete models to judge the accuracy of those models. In addition, the constraint 
can be used as an equation (set) to improve the geopotential. For this purpose 
I have used the (q = 0, second order effects) values in Table 2 divided by 
/2 1Q~^ / ^ and normalized to the maximum, and find a lumped coefficient 
equivalent to Equation (20) to be: 

10® (C, 0 = (- 11-2 ± 1-5. - 32.8 ± 1.5) = (C, 8 ) 15 , - .292(C. S)i, 

— — .803(C, S)jg — .803(C, 3 ) 2 ^ 1 ^ - .565(C, 5)23 

- .286(C S> 2 s, 14 - .0S4(C, S>27, 14 + -100(0, (22) 

The solutions of this equation are plotted in Figure 2. A number of remarks 

on this figure are in order: 

1 . The observed value is in fair agreement with values calculated from 
recent comprehensive fields determined without use of the ERTS orbit. 

2. All the fields displayed use surface gravity and satellite tracking data 
with the exception of GEM 1 which uses only satellite optical data 
(Smith, Lerch and Wagner, 1973). This aspect does not seem too 
significant except that the "best" result is achieved for the WG^ 72 
field (Department of Defense, classified, 1974) which combines 
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considerably more data than the others. The GEM 6 model (Lerch, et al., 
1974) contains considerably more satellite data than GEM 1, as well as 
surface gravity information, yet shows almost no improvement in pre- 
dicting the ERTS nodal variation. The gravitational parameters of WGS-72 
are determined from extensive Doppler laser and optical tracking on 
about 30 satellite orbits combined with surface gravity and astrogeodetic 
data. The field is complete to (19, 19) with resonant and zonal terms as 
high as 28th degree. The SAG SE 3 (Gaposchkih, 1973) though it uses 
more satellite (especially laser) and surface data than SAO SE 2 
(Gaposchkin and Lambeck, 1970) gjives only a marginally Improved ERTS 
result. 

3. While taken together the fields show considerable variability (up to 32% 
in amplitude) from the ERTS result, the variation is much smaller 
compared to the amplitude estimated on the basis of Kaula's rule. 

This estimation is /2 times the root sum of the squares of the terms 
in (22) with (C, S)^ = ICT^ . The relevant terms of WGS-72 for 

example are (on average) 74% of the rule, while the predicted amplitude 
with WGS-72 is only 34% of Kaula’s rule. The variability of the 
amplitudes for all the fields is only 14% with respect to the rule. 

Finally, it should be emphasized again that these second order effects, 
proportional to (n//r)^ , are particularly Important for the variations in deep 
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resonaBce (commensur ability) where 'Z' “■ 0. For example, they do not appear to 
have been considered by King Hele (1974) in his analysis of the node and perigee 
of the orbit of Cosmos 72 during its resonance pass with 15th order geopotential 
terms in 1972. The character of the perturbation through commensurability 
for the node and perigee remains to be investigated. These variables probably 
do not librate since the acceleration of is under the control of the stronger 
resonance in the mean motion. 

SUMMARY 

The interaction of secular oblateness effects with long periodic resonance 
on satellite orbits is calculated. The result is to add a perturbation (of the 
same phase and period) of order ^ (n/^) to the direct (first order) 

resonance effect on the node, perigee and mean anomaly. The direct effect is of 
order {n/p ). Where 0 is small (shallow and deep resonance) the interaction 
term is significant or dominates the resonant motion of the node and perigee. 
These effects have not been considered in previous analytic orbit analyses of 
geodetic satellites. For a number of shallow resonant earth orbits this failure 
results in significant error with respect to the cross track information from the 
node. The effects are especially large on the node of the ERTS orbit (18 day 
period) for viduch the first order theory wovild underestimate the amplitude (30m) 
by a factor of 2. 

Results from the ERTS orbit yield a geopotential constraint which will be 
useful in refining current models for 14th order terms to about the 29th degree. 
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APPENDIX 

INCLINATION FUNCTION RATIOS 

The dominance of direct or indirect terms in the resonance of the node, 
perigee and mean anomaly depends critically on the ratio (3 F/3 1)/F of the 
derivative of the inclination function to its value- For medium to high inclina- 
tion satellites the resonance is generally dominated by the terms of a species 
where - m is small and q = 0 since they are of the lowest degree and do not 
contain the eccentricity as a factor (Allan, 1965; Allan, 1973). 

It is instructive (and simplest) to work out the root mean square (rms) value 
of this ratio for these dominant functions where = m, q = 0 and k = 1; the odd 
order resonances. Allan (1973) gives the normalized inclination function for 
these as: 


F = 




2“ |-(m- 1) ! |(m+ 1) 


cm+lsH>-l 


where 


=2(2m+l)/(2m)!. 


and 


c = cos 1/2, s = sin 1/2 


The full inclination function for the direct term in (14) is: 


F' = (^F/31)/F 
sin I 


(Al) 
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Expanded, tMs is simply: 


F' = 1) + i (m - 1) s“’' 2 cm+ 2 '| / 

sin I L 2 2 J/ 




= _J™ |(m + 1) c“^s + ^ (m - 1) cs"^ = . - 4 (n' + 1) 

sin iD ^ J sinIcsL2 


+ - (m - 1) c2 


y _ 


2 s in I cs 


[m(c^ - s^) - (c^ + s^)] 


= [jn cos I - 1] /s in^ I 


Except for m = 1, F’ is singular at I = 0, 180®. But since both (BF/BI)/sin I and 
F are zero at I = 0, 180® for m > 2, both direct and indirect effects are negligible 
for the close circular orbit equatorial satellite. 

Therefore only the range of 30®- 150° will be taken for the rms computation. 
The indefinite integral: 


3 


cot^ I 1 

3 sin^ I ' 3 sini 


- cos I 2 cot I 
3 s i n^ I ^ 


Therefore: 


aS7t/6 

rms F' (30® - 150°) = A (F')2dl 

2tt j , 

•'Tr/6 


F3 

= i (3.46m + 6.93) 
2tt j 


(A2) 
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For the 13 rev/day orbit, Equation (A2) yields: 

rms - 150°) = 4.98. 

This compares to the value 4 estimated in the text from a cruder sinusoidal 
model for F. 
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Table 1 

Resonance Perturbations of the Node with Second Order Effects 
Effects Due to = /2 10 ^ 

'Urn 


Satellite 

Orbit 

Inclination 

<deg's) 

Apogee 

Ht. 

(km) 

Eccentricity 

Maximum Resonance 
Perturbation Term; meters 

Resonant 

First Order 
Only 

With Second 
Order 

Period 

(days) 

Order 

Agena (1964-lA) 

69.9 

926 



10.8 

5.0 

14 

Anna IB <1962 -60A) 

50.1 

1184 



5.3 

4.8 

13 

BEB (1964-64A) 

79.7 

1075 

.0135 


5.4 

3.0 

14 

BE-C (1965-32A) 

41.2 

1322 

.0257 


4.8 

5.6 

13 

Courtier IB (1960-13A) 

28.3 

1211 

.0161 

1.4 

1.4 

3.8 

13 

DI-C (1967-llA) 

40.0 

1354 

.0532 

2.8 

2.9 

2.5 

14 

DI-D (1967-14A) 

39.5 

1890 

.0848 

8.7 

9.2 

8.4 

13 

Echo Rocket (1960-9B) 

47.21 

1682 

.0118 

65.3 

59.6 

11.9 

12 

GEOS 1 (1965-89A) 

59.4 

2277 

.0719 

6.2 

7.2 


12 

GEOS 2 (1968-2A) 

105.8 

1588 

.033 

5.9 

4.4 

5.7 

13 

GRS (1963-26A) 

49.8 

1294 

.0598 

20.6 

25.3 


14 

Injun 

66.8 

995 

.0079 

7.8 

7.6 

3.8 

14 

Midas 

95.8 

3728 

.0112 

0.9 

0.8 

3.0 

9 

OGO 2 

87.4 

1515 

.0752 

14.0 

14.3 

3.8 

14 

Oscar '7 

89.7 

1199 

.0224 

4.7 

4.9 

2.2 

14 

Secor 5 

69.2 

2420 

.0793 

3.2 

3.2 

3.4 

12 

Transit 

66.8 


.0076 

7.2 

7.0 

3.5 

14 

Vanguard 2 

32.9 

3282 

.1641 

3.3 

3.2 

2.7 

11 

8BN-2 

90.0 

1128 

.0058 

2.8 

3.1 

2.4 

13 

Vanguard 3 

33.3 

3748 

.1901 

24,9 

42.2 

12.2 

11 

Telstar 

44.8 

5640 

.2429 

22.9 

21.8 

14.9 

9 

EBTS-1 

99.1 

908 

.0080 

10.2 

40.2 

18.0 

14 
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Table 2 

Resonant Effects on the Node of the ERTS Orbit 
Unit = 10”'* Degrees 
Amplitude of Effects Due to (C, S)^ m ~ 


('i, m) 
q = 0 

First Order 
Effect Only 

With Second 
Order Effect 

(■1, m) 

l\dth Second 
Order Effect 

(^,m) 
q = +1 

With Second 
Order Effect 

15,14 

-0.72 

-3.71 

14,14 

.00 

14,14 

-.02 

17,14 

1.91 

.85 

16,14 

.02 

16,14 

.03 

19,14 

2.05 

1.86 

18,14 

.00 

18,14 

.03 

21,14 

1.41 

1.52 

20,14 

-.01 

20,14 

.02 

23,14 

.72 

.90 





25,14 

.24 

.38 





27,14 

-.03 

.06 





29,14 

-.15 

-.10 





31,14 

-.18 

-.15 





28,28 

.05 






30,28 

.05 






32,28 

CO 

o 

* 

1 






34,28 

-.05 






36,28 

-.04 
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O OBSERVED MEAN ELEMENT 
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FROM A TRAJECTORY WITHOUT 
RESONANCE 


ANALYTIC VARIATION 

WITH SECOND ORDER EFFECTS 

WITH ONLY FIRST ORDER EFFECTS 



TIME: MODIFIED JULIAN DATE (DAYS) 


Figure 1. Resonant variation of the node of the ERTS orbit. 





Figure 2. Lumped coefficient from resonance in node for ERTS. 



